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O . ABSTRACT 

■ We build template maps for the polarized Galactic-synchrotron emission on large 

angular scales (FWHM = 7°), in the 20-90 GHz microwave range, by using WMAP 
data. The method, presented in a recent work, requires a synchrotron total intensity 
survey and the polarization horizon to model the polarized intensity and a starlight 
polarization map to model polarization angles. The basic template is obtained directly 
at 23 GHz with about 94% sky-coverage by using the synchrotron map released by 
the WMAP team. Extrapolations to 32, 60 and 90 GHz are performed by computing 

Q . a synchrotron spectral index map, which strongly reduces previous uncertainties in 

^ ■ passing from low (1.4 GHz) to microwave frequencies. Differing from low frequency 

, data, none of our templates presents relevant structures out of the Galactic Plane. Our 

■ map at 90 GHz suggests that the synchrotron emission at high Galactic latitudes is 

low enough to allow a robust detection of the £^-mode component of the cosmological 
signal on large-scale, even in models with low-reionization (r = 0.05). Detection of 
the weaker B-m.ode on the largest scales < 10) might be jeopardized unless the 

I value r = 0.17 found by WMAP is confirmed, and T/ S > 0.1. For lower levels of the 

gravitational-wave background the B-mode seems to be accessible only at the ^ ~ 100 
peak and in selected low-synchrotron emission areas. 

Key words: polarization. Galaxy, cosmic microwave background. Method: numerical. 



1 INTRODUCTION 

Recent results by DASI (Kovac et al. 2002) and WMAP 
(Bennett et al 2003a) have opened a new window in studying 
both the Cosmic Microwave Background (CMB) and the 
Galaxy. 

The polarized signal detected with the DASI ex- 
periment, as well as the WMAP measurement of the 
temperature-polarization cross-spectrum Cf^, emphasized 
the importance of studying CMB Polarization (CMBP). In 
particular, the discovery of a relevant signal on large angu- 
lar scales performed by WMAP has provided evidence of an 
unexpectedly early reionization era (Kogut et al. 2003). 

Furthermore, WMAP provided the first all-sky total- 
intensity maps of the microwave Galactic emission, allowing 
an insight into synchrotron, dust and free- free components 
of the Galaxy (Bennett et al. 2003b). 



A sound measurement of CMBP requires good knowl- 
edge of polarized foregrounds (both Galactic and extra- 
galactic), which are potentially more dangerous than in to- 
tal intensity. This is even more true for the B-mode, re- 
taining information on the gravitational-wave background 
(Kamionkowski & Kosowski, 1998), whose emission level is 
expected to be orders of magnitude below that of the E- 
mode. 

Among all foreground components, the Galactic syn- 
chrotron radiation is expected to be the most relevant up 
to 100 GHz. The lack of data in both the CMB frequency 
range and at high Galactic latitudes makes the building of 
templates necessary (Kogut & Hinshaw 2000, Giardino et 
al. 2002, Bernardi et al. 2003a, hereafter B03). These can 
substantially help in studying how to extract CMB maps 
from contaminated signals (e.g. see Bennett et al. 2003b). 
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In BOS we presented a method modelling the Galac- 
tic polarized-synchrotron emission by using the radio- 
continuum total intensity surveys available at low frequency 
(< 1.4 GHz) to model the polarized intensity, and starlight 
polarization optical data as a template for polarization an- 
gles. 

Template maps obtained with this method are virtually 
free of Faraday rotation. However, when extrapolating low 
frequency data (the only available before the latest WMAP 
release) to the cosmological window, uncertainties arise on 
the spectral index to be used. In BOS we used the mean spec- 
tral index derived by Platania et al. (1998) up to 19 GHz, 
since no information were available at higher frequencies. In 
addition, variations across the sky were too poorly known 
to be taken into account. 

The WMAP release of total-intensity synchrotron maps 
in the cosmological window gives us the possibility to signi- 
ficatively improve our templates. Actually, the application 
of our method directly at frequencies of interest for CMBP 
fully avoids extrapolation uncertainties. 

Our templates are built with FWHM = 7° at the three 
frequencies of 2S, S2 and 90 GHz interesting for the SPOrt 
experiment (Cortiglioni et al. 2004), one of whose aims is 
studying the CMBP foregrounds. We also generate a tem- 
plate at 60 GHz, this frequency being of interest for other 
CMB experiments (WMAP, Planck). Ah of the templates 
cover large fractions of the sky, namely 94% at 2S GHz 
and 92% at S2, 60 and 90 GHz. In addition, we calcu- 
late polarized-syncrotron angular power spectra both for 
the full-sky templates and for the high Galactic latitude 
emission, discussing possible effects on the measurement of 
CMBP E and 5-mode spectra. 

The paper is organized as follows: in Section 2 we apply 
our method to WMAP data at 2S GHz and provide a Galac- 
tic polarized synchrotron template at the same frequency; 
in Section S we extrapolate our results up to 90 GHz; in 
Section 4 we describe the angular power spectra of our full- 
sky templates at the different frequencies, in Section 5 we 
compute the angular power spectra for the high Galactic 
latitude emission and compare synchrotron to CMBP and, 
finally, conclusions are presented in Section 6. 



2 A POLARIZED TEMPLATE AT 23 GHZ 
USING WMAP DATA 

In BOS we presented a method to obtain a template for the 
polarized Galactic synchrotron emission based on total in- 
tensity radio maps, starlight polarization data and the mod- 
ellization of a polarization horizon as suggested by several 
authors (Duncan et al. 1997, Gaensler et al. 2001, Landecker 
et al. 2002). 

The method is based on separate derivations of the po- 
larized intensity and polarization angle maps. Whereas the 
latter is obtained from starlight polarization data, the for- 
mer is the result of the application to total-intensity syn- 
chrotron maps of a filter mimicking the effects of the polar- 
ization horizon. 

The main steps of our method can be summarized as 
follows: 

(i) a total intensity synchrotron map is obtained from 



available surveys, after separating the CMB and free-free 
components when necessary; 

(ii) a filter accounting for the effects of the polarization 
horizon is appli ed to the / map, providing the polarized 
intensity Ip — yQ^^+lP] 

(iii) a polarization-angle map is obtained from starlight 
data by interpolating the non-uniform data coverage; 

(iv) the polarized-intensity and polarization-angle maps 
are used to compute Q and U templates including a con- 
volution with a FWHM — 7° Gaussian filter to match the 
SPOrt angular resolution (finer resolutions are not possible 
due to the sampling of starlight polarization data). 

Before WMAP, total intensity surveys were available up 
to 1.4 GHz and our templates were limited by our ignorance 
of the spectral index distribution needed to extrapolate to 
higher frequencies. Indeed, previous templates at 22, S2, 60 
and 90 GHz were obtained by using a constant spectral index 
(3 = — S, corresponding to the mean value between 1.4 and 
19 GHz (Platania et al. 1998). 

The model is now substantially improved by using the 
latest WMAP data, including a total-intensity synchrotron 
map (cleaned from CMB, free-free and dust contributions) 
with resolution of 1° at 2S GHz, that directly allows us to 
compute the template at a frequency interesting for CMBP 
purposes. 

The polarized intensity map is obtained by applying to 
the 23 GHz synchrotron map the equation 

i,('^,i,b)=p^i{i^,i,b) (1) 

mimicking the polarization horizon effects (see BOS). Here 
Ip and / are the polarized and total intensity of the syn- 
chrotron radiation, respectively, L{l,b) is the thickness of 
the synchrotron emitting region, p is the polarization de- 
gree, Rph is the distance of the polarization horizon. We 
adopt the same normalization factor pRph = 0.9 Kpc as 
computed in BOS. 

From this polarized brightness temperature, and using 
the same polarization-angle map as in BOS, Q, U and Ip 
maps are obtained. They are shown in Figure 1, where the 
grey region corresponding to the North Celestial Pole is an 
area where the interpolation of starlight polarization angles 
failed (see BOS for details). It is worth emphasizing that the 
new maps cover almost 94% of the sky, to be compared with 
the half-sky coverage of the BOS template. 

Confirming the BOS results, the distributon of the po- 
larized intensity is quite different from that of total intensity 
(Figure 2). The latter peaks towards the Galactic centre and 
traces well the Galactic Plane. Instead, the polarized emis- 
sion has a sub-dominant Galactic Centre and the Galactic 
Plane shows a patchy structure (Figure S). 

A comparison with the BOS template shows a general 
agreement between the two maps along the Galactic Plane, 
and in particular in the following three regions: 

(i) Fan region located at 100° < I < 160°; 

(ii) the region at S0° < / < 45°; 

(iii) the region around / ^ 15°, though only partially 
present in the BOS template. 

A comparison with the Parkes survey (Figure 4) shows re- 
markable similarities as well. The brightest polarized region 
in our template is at Galactic longitude S15° < I < S45°, 
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Figure 1. Our 23 GHz polarized synchrotron template obtained using WMAP data as total intensity emission. Maps are convolved with 
a FWHM = 7° Gaussian filter. 
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Figure 2. WMAP's total intensity synchrotron emission at 
23 GHz (from Bennett et al. 2003b) smoothed with a 7° beam. 



where also Duncan et al. (1997) found a high level in the 
polarized background at 2.4 GHz (see also Figure 4). The 
bright features present in our template correspond to emis- 
sion spots in the Duncan map, with the exception of the 
large region centred at / ^ 300°. This discrepancy is the 
same noted in B03 and is not surprising. Actually, this is a 
peculiar area with strong variations in Rotation Measures 
(Sofue & Fujimoto 1983, Duncan et al. 1997, B03), so that 
the absence of polarized emission in the Parkes data is likely 
to be due to Faraday depolarization. 

Out of the Galactic Plane, the model computed with 
the WMAP data (hereafter the 23-GHz-W template) does 
not show any relevant structure. This differs from the emis- 
sion at 1.4 GHz, where the Northern Galactic Spur is among 
the main polarized structures, both in the Brouw & Spoel- 
stra (1976) data and our 1.4 GHz B03 model. The situation 
is similar in the total-intensity case, where the Northern 
Galactic Spur is a prominent feature at 1.4 GHz, but is no 
more significant at 23 GHz. This behaviour could be ex- 
plained following Bennett et al. (2003b), who showed that 
the spectral index of the synchrotron radiation in the 0.408- 
23 GHz range is quite different in the Plane and out of it. 
The Galactic Plane having a flatter spectral index becomes 
more and more dominant with increasing frequency. 



3 SPECTRAL INDEX BEHAVIOUR AND 
EXTRAPOLATION TO HIGHER 
FREQUENCIES 

The most important frequency range for CMBP studies is 
around 90 GHz, where the polarized foregrounds are ex- 
pected to have minimum emission. Since the 94 GHz syn- 
chrotron map provided by WMAP appears to be dominated 
by noise at high Galactic latitudes, we do not apply our 
method directly to it to provide a synchrotron template at 
90 GHz. Instead, we prefer to extrapolate our 23 GHz tem- 
plate, characterised by a better S/N ratio. 

Faraday rotation effects are negligible at frequencies 
higher than 20 GHz (e.g. see B03). In absence of RM modu- 
lation, polarization and total intensity are expected to have 
the same behaviour. 

Synchrotron spectral-index variations across the sky 
have to be carefully modelled when extrapolating syncrotron 
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Figure 3. Comparison between the total (observed, solid line) 
and polarized (template, dashed line) brightness temperatures 
along the Galactic Plane at 23 GHz. We consider all the pixels 
with |6| < 0.1° 



I (Parkes survey) 




Figure 4. Polarized intensity emission from the Parkes survey at 
2.4 GHz smoothed on 7° (Duncan et al. 1997). 



maps to higher frequencies. Bennett et al. (2003b) provided 
a spectral-index map, for the 0.408-23 GHz range, which 
applies to a mixture of synchrotron and free-free emissions 
and is not directly useful for our purposes. 

We thus extract a nearly full-sky synchrotron map at 
1.4 GHz, by means of a modified Dodelson technique (B03) 
to separate between components, from the Haslam et al. 
(1982) survey at 0.408 GHz, the Reich (1982) survey at 
1.4 GHz, and the Jonas et al. (Jonas, Baart & Nicolson, 
1998) survey at 2.3 GHz. 

Properties of data at 0.408 GHz and 1.4 GHz have al- 
ready been described in B03 (and references therein). Data 
at 2.3 GHz correspond to a radio-continuum survey of the 
Southern sky (—83° < ^ < 13°), with angular resolution of 
20 arcmin. 

All maps are smoothed to the same resolution of 51 ar- 
cmin, corresponding to data at 0.408 GHz. The free-free 
subtraction is performed adopting the B03 model for the 
spatial distribution of the synchrotron spectral index which, 
being symmetric with respect to both the Galactic Plane 
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Figure 5. Total intensity synchrotron map at 1.4 GHz on 7° 
(cleaned from free-free contribution) obtained from Haslam et al. 
(1982), Reich (1982) and Jonas et al. (1998) maps. 




Figure 6. Spectral index map of the synchrotron emission on 7° 
angular scales between 1.4 and 23 GHz. The grey region around 
the South Celestial Pole is the unobserved area in the 2.3 GHz 
survey. 



and the Galactic centre, can be extended to the Southern 
sky. 

Due to different sky coverages of the available surveys, 
synchrotron maps are built separately in the Northern and 
in the Southern emispheres. Where the two maps overlap, 
we retain information from the Northern one, since the zero 
level of its parent surveys are more accurately determined 
than that at 2.3 GHz. Moreover, we correct the offset of 
the Southern synchrotron map shifting its mean level to the 
that of the Northern map in the overlapping region. Discon- 
tinuities at borders are at negligible levels for our purposes, 
being cancelled out by the convolution with a FWHM = 7° 
Gaussian filter. The final map, covering almost 99% of the 
sky, is shown in Figure 5. 

A comparison with the 23 GHz synchrotron map pro- 
vided by WMAP allows us to determine the synchrotron 
spectral index in the 1.4-23 GHz range (Figure 6). 

The overall pattern of our synchrotron spectral-index 
map is substantially in agreement with that by Bennett et al. 
(2003b): the spectral index is rather fiat along the Galactic 
Plane, whereas steeper spectra are found at high Galactic 
latitudes. As shown in Fig. 7, the typical spectral index is 
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Figure 7. Synchrotron spectral index distribution between 1.4 
and 23 GHz. 



P = —3.1, which is steeper than the value found by Bennett 
et al. (2003b); this might be due to the thermal contribution 
in their data. 

We use our spectral-index map to extrapolate the 
1.4 GHz B03 template up to 23 GHz (hereafter the 23- 
GHz-LF [low frequency] template). The result is shown in 
Figure 8, and can be compared with the 23-GHz-W map to 
further test our template-building procedure. 

A general agreement in the morphology can be ob- 
served. In particular, the most relevant characteristic in 
both maps is the lack of features at high Galactic latitudes. 
The Northern Galactic Spur, present at 1.4 GHz, disappears 
when extrapolated with the spectral-index map of Figure 6, 
in agreement with 23-GHz-W. Also the main features in the 
Galactic Plane are similar. 

Indeed, the model applied to the 0.408-1.4 GHz data 
(including the free-free component separation) and extrap- 
olated to 23 GHz by using spatially varying spectral indeces 
presents results very similar to the 23-GHz-W template, 
which is free from extrapolation uncertainties. 

In order to produce template maps up to 90 GHz, the 
spectral index map in Figure 6 is not sufficient. A steepening 
at higher frequencies, in fact, has been found by Bennett et 
al. (2003b), who reported a mean steepening of 0.5 in the 
23-41 GHz range. To extrapolate our 23-GHz-W template, 
thus, we increase the spectral indexes by this amount. 

Figure 9 shows /p, Q and U maps at 90 GHz, which 
is the most interesting frequency for CMB measurements. 
Maps at 32 and 60 GHz look qualitatively rather similar. 
As already pointed out in this section, the flatter Galactic 
Plane component becomes more and more dominant with 
increasing frequency. A further important characteristic is 
that the Galactic Plane essentially retains its morphology. 
However, again due to the different spectral behaviour, the 
Fan region is fainter and fainter with increasing frequency 
respect with the 300° < / < 330° region. 



4 ANGULAR PROPERTIES OF THE 
POLARIZED TEMPLATE 

The distribution of the polarized emission at different angu- 
lar scales is described by the angular power spectra (APS) 
Cf , Cf and Cf = Cf + Cf . We compute them by the 
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Figure 8. Maps of our template obtained from the 1.4 GHz B03 one scaled to 23 GHz (23-GHz-LF). The angular resolution is 7°. 
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Figure 9. Maps of our template at 90 GHz. The angular resolution is 7°. 
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Figure 10. Cf^ Cf, Cf and C^^ spectra computed for the 23-GHz-LF (solid) and 23-GHz-W templates (dotted). Spectra for the 
32 (dashed), 60 (dashed-dotted) and 90 GHz (dashed-triple dotted) models are shown as well. 



Table 1. Best fit APS slopes and amplitudes as obtained from our templates in the 2 < £ < 25 range using WMAP data. First and 
second lines report the 23-GHz-LF and 23-GHz-W templates respectively. 



u as C^Q as ap C[q aj^ C^g 

(GHz) a^K2) 0^K2) (/iK^) 

23-LF 1.34 ±0.09 37 ib 2 1.57 ±0.09 58 ± 4 1.54 ± 0.09 101 ± 6 2.00 ± 0.09 115 ± 7 

23-W 1.55 ±0.09 64 ±4 0.83 ± 0.09 49 ± 3 1.48 ± 0.09 148 ± 9 1.00 ± 0.09 123 ± 8 

32 1.54 ±0.09 9.7 ±0.6 0.79 ± 0.09 7.3 ± 0.4 1.43 ± 0.09 22 ± 1 1.00 ± 0.09 20 ± 1 

60 1.47 ±0.09 0.27 ±0.02 0.65 ± 0.09 0.20 ± 0.01 1.31 ± 0.09 0.60 ± 0.04 0.94 ± 0.09 0.60 ± 0.04 

90 1.41 ±0.09 (27 ± 2) X 10-3 0.58 ± 0.09 (20 ± 1) x 10"^ 1.23 ± 0.09 (60 ± 4) x 10"^ 0.90 ± 0.09 (65 ± 4) x 10" 



correlation-function method (see Sbarra et al. 2003 for de- 
tails), accounting for irregularities in the sky coverage. The 
APS of our templates are reported in Figure 10. Their over- 
all behaviour, with particular reference to , can be rep- 
resented by a power law: 

Cj = C7jo (^) Y = E,B, P, 7p, (2) 

whose best fit parameters in the 2 < ^ < 25 range are re- 
ported in Table 1. We perform a linear fit to the quantities 
log d and log i giving greater weights to high order multi- 
poles (since Ci is an average of (2i -\- 1) squared harmonic 
amplitudes). In addition, we compute the scalar spectrum 
C/ of Ip for comparison with other works. 



The two templates at 23 GHz are in good agreement in 
the overall intensity, whereas the slopes are somewhat dif- 
ferent, the 23-GHz-W being flatter even if the slopes of 
the total polarization spectrum are statistically equal. 
Differencies might be due to the different sky coverages of 
the two maps. In fact, we have computed the APS for the 
23-GHz-W template hmited to the coverage of 23-GHz-LF 
and have found a better agreement between the two tem- 
plates. This suggests the Southern sky contributes to make 
the overall slope flatter. 

A slight decrement of the slopes with increasing fre- 
quencies appeares in the Cf and Cf spectra of our tem- 
plates, the 90 GHz ones being the flattest. However, these 
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variations are within 2a of the slopes themselves and a ro- 
bust result cannot be claimed. 

While the slopes of Cf and are compatible with 
each other, Cf and C^^ are significantly flatter and this 
behaviour is still not clear. In addition, Cf and C[ slopes 
are in agreement with those previously calculated in real 
radio surveys (Tucci et al. 2000, Baccigalupi et al. 2001, 
Bruscoli et al. 2002, Giardino et al. 2002). However, one 
should take into account that a perfect agreement is not 
expected at all, since our template and previous works differ 
in both the angular scale and the frequency (1.4-2.7 GHz 
for the real surveys). 



5 HIGH GALACTIC LATITUDE EMISSION 

The templates presented here predict that the synchrotron 
polarized emission is concentrated on the Galactic Plane. As 
it already occurs for the total intensity, this suggests us to 
search for the CMBP signal at high Galactic latitude. 

Let's us recall that this was not that obvious before 
the present work, since large and bright features exist, out 
of the Galactic Plane, in low frequency observations (see, 
for example, the Northern Galactic Spur in the 1.4 GHz 
polarized data of Brouw & Spoelstra 1976). 

To study the expected synchrotron emission at high 
Galactic latitudes we mask off the Galactic Plane by using 
the pixel-count histogram of the Ip map at 23 GHz shown 
in Figure 11, in a way similar to that described by Bennett 
et al. (2003b). The intensity distribution is strongly asym- 
metric, and pixels with the higher values correspond to the 
Galactic Plane emission to be cut away. To set the value 
of the threshold defining which pixels are to be discarded, 
we fit the peak of the distribution and subtract the best-fit 
curve from the histogram. The peak of the residual distribu- 
tion was chosen by Bennett et al. (2003b) as the threshold 
value. 

Differing from the Gaussian behaviour followed by tem- 
perature data (Bennett et al. 2003b), the low emission part 
of the polarized intensity Ip = yQ^^+lP' is better repre- 
sented by a Rayleigh distribution, as it would be if Q and 
U were gaussian distributed: 

F{Ip)^AIpe-'v/^^^") (3) 

where A and a are parameters to be determined. 

We fit the distribution in the 0-22 /xK range. The best- 
fit subtracted histogram has a peak at 39.1 /iK; however, 
we prefer to be more conservative and set the threshold at 
34.1 /iK of the secondary peak. The pixel mask of the re- 
tained pixels {Ip < 34.1 /xK) covers about 68% of the sky 
and is shown in Figure 12. 

APS computed on our template maps after applying 
this mask are shown in Figure 13. Power spectra at high 
Galactic latitudes can be represented by power laws apart 
from the region at large where they are observed to flatten, 
probably due to pixel noise. Therefore, we fit the 2 < ^ < 15 
range and report the results in Table 2. 

At each frequency, power spectra levels are about two 
orders of magnitude below those including the Galactic 
Plane which, therefore, dominates the APS calculated in 
Section 4. 

Slopes are steeper than those computed on the overall 
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Figure 11. Histogram of the polarized intensity distribution in 
the 23-GHz-W template. The solid curve represents the Rayleigh 
function used to fit the distribution e the dashed line represents 
the residual distribution. 
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Figure 12. The pixel mask in Mollweide projection. Grey pixels 
represent the excluded region. 



templates, in particular for Cf , suggesting that the flat- 
tening observed in the full-sky templates is mainly due to 
Galactic Plane features. No significative flattening with in- 
creasing frequency is found at high Galactic latitude. 

Figure 14 shows a comparison between our estimate of 
synchrotron power spectra at high Galactic latitudes and 
expected CMBP angular spectra. The power spectra of our 
templates are converted to CMB thermodynamic tempera- 
ture by multiplying them by the square of the conversion 
factor 

where x = hu/kTcmb ~ z^/56.8 GHz. CMBP angular spectra 
are computed assuming cosmological-parameter values as in 
the concordance model determined by WMAP data (Spergel 
et al. 2003): ho 71, 0.044, Qcdm 0.226, Qa 

0.73, T = 0.17 and initial adiabatic fluctuations. 

Our templates suggest that at 90 GHz CMBP ^-mode 
measurements in the 2 < ^ < 20 multipole range should 
be only marginally contaminated by synchrotron polarized 
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Table 2. Best fit APS slopes and amplitudes as obtained from our templates in the 2 < £ < 15 range after the pixel mask is applied. 



(GHz) 



^10 



OLB 



^10 



dp 



^10 



^10 



23 1.8 ±0.2 1.3 ±0.1 1.6 ±0.2 

32 1.9 ±0.2 0.13 ±0.01 1.5 ± 0.2 

60 1.9 ± 0.2 (1.7 ± 0.1) X 10-3 1.4 ± 0.2 

90 1.9 ±0.2 (1.04 ± 0.07) X 10-^ 1.3 ± 0.2 



1.4 ±0.1 1.7 ±0.2 

0.13 ±0.01 1.6 ±0.2 

(1.7 ±0.1) X 10-3 1.5 ±0.2 

(1.04 ±0.07) X 10-^ 1.5 ±0.2 



2.6 ±0.2 
0.25 ± 0.02 
(3.3 ±0.2) X 10-3 
(2.0 ±0.1) X 10-^ 



1.5 ±0.2 
1.2 ±0.2 
1.2 ±0.2 
1.2 ±0.2 



1.2 ±0.1 
0.13 ±0.02 
(2.1 ±0.1) X 10-3 
(1.4 ±0.1) X 10-4 
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Figure 13. Cf, , and C^^ spectra computed for the 23 GHz template with the high Galactic latitude mask applied (dotted 
line). Spectra for the 32 (dashed), 60 (dashed-dotted) and 90 GHz (dashed-triple dotted) models are shown as well. 



radiation. Furthermore, even a less reionized Universe (r = 
0.05) should be accessible at this frequency. 

It is worth noticing that even relaxing our hypotheses 
of a steepening by 0.5 in the spectral indeces between 23 
and 90 GHz, e.g. simply adopting the spectral-index map of 
Figure 6, the normalization of the synchrotron power spec- 
trum would be just a factor ~ 4 higher, which does not 
substantially modify the conclusions drawn above. In fact, 
the CMBP power spectrum of models with r = 0.17 and 
r — 0.05 are about two and one orders of magnitude larger 
than the synchrotron spectrum, respectively. 

Obtaining a clean measurement of the CMBP B-mode 
looks more difficult. Large values of the optical depth r push 
a large-scale peak at ^ < 10 and open a second window to de- 
tect this faint signal beyond that at £ ~ 100 (Kamionkowski 
& Kosowsky 1998). Figure 14 shows two models with the 



same cosmological parameters as above, but with the ad- 
dition of a gravitational wave background with tensor to 
scalar fluctuation power ratio T /S = 0.1 and T/S = 0.01, 
respectively. While the former looks reasonably free from 
synchrotron contamination, at least in the low multipole 
range {£ < 10), the level of the latter is comparable with 
our estimate of the synchrotron spectrum. 

We can thus conclude that, for values r = 0.17 and 
T/S > 0.1, a CMBP B-mode measurement on large angular 
scale free from synchrotron contamination is achievable. For 
lower T/S* values, instead, detection of the B-mode would be 
only possible in selected regions with low foreground emis- 
sion, and large enough to access the i ~ 100 peak. A good 
example is the patch observed by the BOOMERanG exper- 
iment and selected for BaR-SPOrt (Cortiglioni et al. 2003), 
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Figure 14. Left: Cf computed on the 23 GHz (dotted), 32 GHz (dashed), 60 GHz (dashed-dotted) and 90 GHz models (dashed-triple 
dotted). The sohd and long-dashed lines represent the CMB ^-mode spectrum of the concordance model with r = 0.17 and r = 0.05, 
respectively. Erros bars represent the cosmic variance. Right: the same but for . The solid and long-dashed lines represent CMB 
S-mode spectrum with T/S = 0.1 and T/S = 0.01, respectively 



Table 3. Prms of our templates at the four reference frequencies 
in thermodynamic temperature. The -Prms is computed on the 
remaining 68% of the sky after applied the pixel mask. Bottom 
line reports the Prms of the CMBP E-mode model with r = 0.17. 
The error on the CMBP Prms takes into account only the cosmic 



contribution. 

1/ (GHz) 




Prms (mK) 


23 




10 


32 




3.2 


60 




0.36 


90 




0.095 


CMBP E-mode (r 
CMBP ^-mode (r 


= 0.17) 
= 0.05) 


0.49 ± 0.03 
0.189 lb 0.007 



for which synchrotron contamination is expected to be under 
control already for T/S > 0.01 (Bernardi et al. 2003b). 

To complete our comparison between synchrotron and 
CMB, we compute the Prms of our masked template maps 
at the four frequencies (Table 3) using the equation (Zaldar- 
riaga 1998) 

Prms = (5) 

where Wi is the window function defined as 

Wi = e-'^'^'^-" (6) 

and cr = FWHM/A/8hr2 

At 90 GHz, the Prms signal is already lower than the 
value computed in B03 using the half faintest sky. Consider- 
ing the 50% faintest sky, a Prms ^ 0.057 /xK is obtained with 
the present model, showing that the mean polarized value of 
the synchrotron background appears quite low in large part 
of the sky. 

Table 3 shows that a r = 0.17 model has a Prms for the 
CMB E-mode which is about five times the Prms of the high 
latitude synchrotron foreground. This suggests that even an 



experiment having just the sensitivity to detect the Prms has 
good chances to measure the optical depth r. 



6 CONCLUSIONS 

In this work we generate template maps of polarized Galac- 
tic synchrotron emission at 23, 32, 60 and 90 GHz using 
WMAP data. 

We apply the method developed in B03 to the 23 GHz 
total intensity synchrotron map released by the WMAP 
team. This greatly improves the B03 templates avoiding the 
uncertainties in the extrapolation from low frequency. As 
for the previous work, the angular resolution is limited to 
FWHM=7° mainly due to the starlight polarization data 
sampling. 

The basic template (23-GHz-W) is constructed at 
23 GHz applying the method directly to the WMAP data 
and it covers almost all the sky (94%). Differing from 
the 1.4 GHz B03 template where the Northern Galactic 
Spur represents an important feature, 23-GHz-W does not 
present relevant structures out of the Galactic Plane. This 
replicates what already occurs in total intensity and leaves 
most part of the sky with low synchrotron emission and use- 
ful for CMBP investigations. 

To achieve the templates at the other frequencies a spec- 
tral index map is required to extrapolate the 23-GHz-W re- 
sult. Differing from Bennett et al. (2003b), which derived a 
spectral index map from 0.408 to 23 GHz with no free-free 
subtraction in the 0.408 GHz map, we computed the fre- 
quency behaviour using pure synchrotron maps at 1.4 and 
23 GHz. The 23 GHz one is that released by the WMAP 
team. The 1.4 GHz map has been achieved from low fre- 
quency surveys by applying a Dodelson-like component sep- 
aration. We use the Haslam et al. 0.408 GHz and the Reich 
1.4 GHz surveys for the Northern sky and the Haslam et al. 
and the Jonas et al. 2.3 GHz ones for the Southern Hemi- 
sphere. The results is the first almost all-sky spectral index 
map of the Galactic synchrotron emission in the 1.4-23 GHz 
range. Furthermore, we add an extra 0.5 to the indeces to 



12 G. Bernardi et al. 



account for the steepening found by Bennett et al. (2003b) 
in the 23-41 GHz range and this spectral index map is used 
to compute template maps at 32, 60 and 90 GHz covering 
92% of the sky. 

The spectral index map provides steeper indexes at 
high Galactic latitudes so that the resulting templates have 
the Galactic Plane emission more and more dominant from 
23 GHz up to 90 GHz. This difference also explains the dis- 
appearing of the Northern Galactic Spur in our templates. 
Unclear, instead, it is the reason for such a flat frequency 
spectrum found in the Galactic Plane, it being either intrin- 
sic to the synchrotron emission (for instance due to a young 
population of relativist ic electrons) or due to a residual ther- 
mal component. 

Beyond the lack of relevant structures at high Galactic 
latitudes, the differences among the templates of this work 
and the B03 one also concern the APS. APS of our full-sky 
templates are dominated by the emission from the Galactic 
Plane and flatter slopes, namely in the 1.4-1.6 range for Cf ^ 
0.5-0.8 for Cf and 1.2-1.5 for Cf . This general flattening 
is still not clear. 

A further unclear point is that the synchrotron APS 
slopes (in particular Cf and Cf) seem to flatten with in- 
creasing frequency. If confirmed, these two aspects would 
confirm that low frequency APS cannot be straightforwardly 
extrapolated from low to microwave frequencies. 

The situation is less complex at high Galactic latitudes. 
In fact, limiting the analysis out of the Galactic Plane, a 
steepening of the slopes occurs, which become more similar 
to the B03 ones, and no significative dependency from the 
frequency is found. 

Interesting conclusions are drawn by the comparison 
of the high Galactic latitude portion of our templates to 
CMBP as from the concordance model of the WMAP first- 
year data. Indeed, the CMBP E-mode spectrum is about 
2 orders of magnitude above the synchrotron signal in our 
template at 90 GHz. Synchrotron is thus unlikely to con- 
taminate the cosmological signal. Furthermore, even a less 
re- ionized Universe (r = 0.05) should be dominant over the 
Galactic emission. 

The situation looks different for the CMBP B-mode. 
Our template predicts that this cosmological signal might 
be accessible at largest scales (£ < 10) only for a Tensor 
to Scalar fiuctuation power ratio T/S > 0.1. To investi- 
gate models with lower values of T/S one should restrict 
to selected low synchrotron emission areas large enough to 
detect the i ~ 100 peak. An example is the sky patch ob- 
served by the BOOMERanG experiment where models with 
T/S > 0.01 seem to be accessible (Bernardi et al. 2003b). 

Our templates obviously rely on the assumption that 
WMAP's dust-correlated signal at 22 GHz is genuine syn- 
chrotron, as supported by Bennet et al. (2003b). However 
several authors (de Oliveira-Costa et al. 2003, Banday et 
al. 2003, Lagache 2003) recently support the opposite view, 
that WMAP's data may be more consistent with anomalous 
dust emission. For instance, while the spinning-dust model 
originally proposed by Draine & Lazarian (1998) may meet 
some difficulty, emission by very small grains is supported 
by Lagache (2003). 

Considering that dust polarization level is likely less 
than 5% rather than the 15-30% of the synchrotron, a mea- 
surement of polarization at 22-32 GHz could help to decide 



on this issue: should the polarization level be found a fac- 
tor ~ 5 below the predictions of our template, the role of 
anomalous dust emission would be supported. 

Also, in such a case, the contamination of CMBP maps 
at 90 GHz would be even smaller than considered in our 
discussion above, and prospects for the detection of a cos- 
mological B mode would be slighly better. 
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